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Abstract. We discuss aspects of open and hidden charm production in deuterium-nucleus collisions at 
RHIC and LHC energies. We describe calculations of the total cc cross section and the charm quark 
transverse momentum distributions. We next explain how shadowing and moderate nuclear absorption 
can explain the PHENIX J/ip dAu/pp ratios and predict the combined effect of shadowing and absorption 
in 6.2 TeV d+Pb collisions. 



1 Introduction 

Heavy flavors were one of the first true hard probes of 
heavy ion collisions through the J/tp measurements of 
NA38 and its successors T'. These same measurements 
also helped show the importance of baseline measurements 
in pp and pA interactions to separate 'normal' from 'anoma- 
lous' behavior. Now, after initial Au-|-Au runs at RHIC, 
the first baseline pp and d-|-Au data have been reported 
for heavy flavors, both open charm and J/iJj. 

Open charm decay leptons were separated for the first 
time in a heavy ion environment by PHENIX in Au+Au 
colhsions at ^/Snn = 130 GeV 0. Now STAR has results 
both in the lepton channel |3] and in the spectra of re- 
constructed D mesons 0]. These measurements will help 
more fully understand open charm production in AA col- 
lisions where energy loss [?], flow jS| and regeneration of 
J/'ip by cc [3 may be important. The J/ipj on the other 
hand, while well measured at the SPS, awaits data from 
the 2004 Au-|-Au run to provide a higher statistics mea- 
surement of J/tp production at 200 GeV. The d-|-Au J/ip 
data, over a wide rapidity range, provide interesting hints 
of the behavior of the nuclear gluon distribution. We dis- 
cuss our dA calculations at RHIC and provide predictions 
for J/tp production in d-|-Pb collisions at the LHC. 



2 Open Charm 

Open charm measurements date back to the late 1970s 
when D and D mesons were first detected, completing the 
picture of the fourth quark begun when the J/ip was de- 
tected in pBe and e"'"e~ interactions. The charm quark was 
postulated to have a mass between 1.2 and 1.8 GeV, calcu- 
lable in perturbative quantum chromodynamics (pQCD). 
Because of its relatively large mass, it is possible to calcu- 
late a total cc cross section, not the case for lighter flavors. 



Charm hadrons are usually detected two ways. The recon- 
struction of charged hadron decays such as I?° — + K^tt'^ 
(3.8%) and D+ K-n+TT+ (9.1%) gives the fufl momen- 
tum of the initial D meson, yielding the best direct mea- 
surement. Charm can also be detected indirectly via semi- 
leptonic decays such as I? — > Klvi although the momen- 
tum of the parent D meson remains unknown. Early mea- 
surements of prompt leptons in beam dump experiments 
assumed that the density of the dump was high enough to 
absorb semi-leptonic decays of non-charm hadrons, leav- 
ing only the charm component. At modern colliders, it is 
not possible to use beam dumps to measure charm from 
leptons but, at sufficiently high px, electrons from charm 
emerge from hadronic cocktails |31IH1- 

Extracting the total charm cross section is a non-trivial 
task. To go from a finite number of measured D mesons in 
a particular decay channel to the total cc cross section one 
must: divide by the branching ratio; correct for the lumi- 
nosity, a£) = No/Lt; extrapolate to full phase space from 
the finite detector acceptance; divide by two to get the 
pair cross section from the single Ds; and multiply by a 
correction factor |U] to account for the unmeasured charm 
hadrons, primarily and Ac- There are assumptions all 
along the way. The most important is the extrapolation to 
full phase space. Before QCD calculations were available, 
the data were extrapolated assuming a power law for the 
xp distribution, related to the longitudinal momentum of 
the charm hadron by xp — Pz/i\^/'2') — 2mT sinh y / ^/S . 
The canonical parameterization for extrapolation over all 
xp is {I — xp)'^ where c was either fit to data over a fi- 
nite xp range or simply assumed. These parameterizations 
could lead to large overestimates of the total cross section 
if < c < 2 was assumed, especially when data were 
taken only near xp — 0. Lepton measurements, partic- 
ularly using beam dumps, resulted in more conservative 
cross sections but were typically at more forward xp. 
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2.1 Total cross section 

There has been a great deal of improvement over the last 
10-15 years. Next-to-leading order (NLO) calculations are 
used in the phase space extrapolation, resulting in con- 
siderably less ambiguity in the shape of the xp distribu- 
tions, da/dxp- The transverse momentum distributions 
are more difficult, as we will discuss later. To calculate 
the total cross section to NLO, scaling functions JOI pro- 
portional to logs of [i} jw? are useful where /i is the scale 
of the hard process. The hadronic cross section in col- 
lisions can be written as 



app{S,m 



E 



PP\~i ■■- I 

dxi dx2 ff{xi,fip) fj{x2,HF) % {s, mI, fJ-n) (1) 

where the sum over i,j runs over q, q and g while xi and 
X2 are the fractional momenta carried by the colliding par- 
tons and ff are the proton parton densities. The partonic 
cross sections are 



a 



/^;^°^(p) + 4^'^'(p)ln(§)]+0(«^)}(2) 



with s the squared partonic center of mass energy, p = 
4to^/s and /fj^'''' are the scaling functions given to NLO 
in Ref. ^Ol- It is most consistent to assume that the fac- 
torization scale, fip, and the renormalization scale, /i^j, 
are equal, fi — jip — Mfl- There is no dependence on the 
kinematic variables. Some NNLO calculations are avail- 
able near threshold, s — X1X2S ~ 1.3 (4m^), applicable 
only for < 20 - 25 GeV [TTIIIS- The NLO correc- 
tions to the leading order (LO) cross sections are relatively 
large, K'-^^ — ctnlo/clo ~ 2 — 3, depending on /i, m and 
the parton densities )13| . 

The NNLO corrections are about as large to next-to- 
next-to-leading logarithm ^Jj but decrease to less than 
K^^^ when subleading logs are included T^. This K factor 
is large because, in the range 1.2 < to < 1.8 GeV, to < 
fj, < 2m with a 5-flavor QCD scale, A^, of 0.153 GeV 
for the GRV98 HO and 0.22 GeV for the MRST parton 
densities, giving 0.21 < as{c) < 0.4. 

Instead of presenting a wide range of possible cross 
sections and emphasizing the uncertainties, the approach 
taken in Ref. [T^ has been to "fit" m and /i for a par- 
ticular parton density and extrapolate to higher energies. 
The results are compared to some of the total cross sec- 
tion data (SI on the left-hand side of Fig. ^ The data 
tend to favor lower values of to, 1.2 — 1.3 GeV. The two 
curves cross each other because the MRST calculation 
with fj, = 2m increases faster at large \/S and smaller x 
due to the stronger QCD evolution of the parton densities 
at the higher scale. Although the fixed target results are in 
good agreement with the calculations, the PHENIX point 
12 at 130 GeV, from Au-f Au electron measurements, and 
the STAR point TS", from a combination of electron and 



reconstructed D measurements, are generally above the 
calculations. The STAR point is about a factor of four 
over the calculated cross section. The higher energy pp 
data from UAl [HI and CDF ^ (not shown, the CDF 
data are, in any case, only at high py, not allowing a total 
cross section measurement) are in better agreement with 
the calculations. 



2.2 Transverse momentum distributions 

Now we turn to the transverse momentum, px, distribu- 
tions. In this case, the quark mass is no longer the only 
scale and pr-dependent logs also appear. Thus, to interpo- 
late between the high pT scale of pr and the low pT scale 
of TO, a scale proportional to toj^, the transverse mass, is 
the natural choice. The charm quark pT distributions are 
not strongly dependent on quark mass for pt > S GeV, 
as may be expected, where the difference in rate is ~ 20% 
between to = 1.2 and 1.8 GeV. The difference in the total 
cross sections is almost all a.t pr < 3 GeV. Changing the 
scale changes the slope of the pT distributions. The distri- 
butions are harder for fi ^ m than fi = 2m. The average 
Pt, (pt), increases with to and is larger for fj, = m. 

It is not enough to equate the charm quark and the 
D meson since hadronization, c ^ D, may involve some 
momentum loss by the quark. If factorization holds in the 
final state (universal fragmentation functions) as well as in 
the initial state (universal parton distributions) then the 
fragmentation functions extracted in e^e~ should also be 
applicable to pp and AA. However, this assumption typ- 
ically does not work well for charm jTH]. The Peterson 
fragmentation function jl9j . typically used in hadropro- 
duction codes, is parameterized as 



D,,h{z) = nUi ^ 



I- z 



(3) 



where Ec = 0.06 was fit to pre-LEP e+e^ data |23, re- 
ducing the charm hadron momentum by 30% relative to 
the charm quark. Current measurements by Belle give a 
slightly smaller value of ec = 0.052 and also suggest that 
this functional form gives the worst fit to the data of all 
the fragmentation functions compared to the data j21| . 
Fragmentation functions calculated in Mellin moment, n, 
space instead of z space, such as the Peterson function, 
tend to predict lower momentum loss by the heavy quark, 
reducing Cc by up to an order of magnitude in the z- 
space representation |22]. (In low collisions, the mo- 
mentum reduction due to fragmentation can be compen- 
sated by intrinsic transverse momentum, kr, broadening. 
However, such broadening cannot compensate the xp dis- 
tributions, only marginally affected by fc-r smearing. We 
have previously shown that the D meson xp distribu- 
tions are consistent with no momentum loss during charm 
quark hadronization The exclusive NLO QQ code 

of Ref. (23) includes Peterson fragmentation and intrinsic 
transverse momentum, fc^, broadening by adding the kx 
kick in the final, rather than the initial state. 



R. Vogt: Open and Hidden Charm Production in dA Collisions at RHIC and LHC 



3 



10' 



10"^ 



10^ 



,4 _ 



10' 



,0 -± 




50 100 500 1000 

Vs" [GeV] 



500010000 



10-= 



10^ 



> 

O 

Xi 10° 

1 

T) 10~l 
\ 

b 

T3 



10 



-2 _ 



10 



-3 





1 

1 




1 1 





/ 


i-V , 








-/ 






















■ w 
1. 








- 


1 


r "\ 




- 


— 




~ " r 























■ ■ V 


















































- - Ix"** 










'— ■ ' >^ 










1 








































' ' ' '^V^ ^ 

■ '" "■ \v 










' - ■> " " 






1 




" "l . , ■ ■ -. 

1 , i- 


I 



10 



Pt (GeV) 



NLO total cc cross sections as a function of x/S (left-hand side) and charm quark pr distribution at \/S = 200 GeV 
\y\ 5: 1 (right-hand side) in pp interactions. The curves are MRST HO (solid) with m — 1.2 GeV and /i^ — 4m^ and 
(dashed) with m = 1.3 GeV and fi^ = m^. The two histograms include fragmentation with the Peterson function 



Fig. 1. The 

in the range 
GRV98 HO 

using Ec ~ 0.06 (dot-dashed) and 0.006 (dotted) 



The effects of fragmentation and intrinsic kx broaden- 
ing of (k^) = 1 GeV^ compensate each other at = 20 
GeV to give a D meson px distribution very similar to 
that of the charm quark. This outcome is desirable be- 
cause the D Pt and xp distributions are similar to those 
of the charm quark at fixed-target energies ^Hl- However, 
at RHIC energies, due to the higher (px) at larger VS, the 
effect of broadening is relatively small and cannot compen- 
sate for the momentum loss induced by fragmentation. In- 
terestingly enough, the STAR D and D* pT distribution 
also agrees rather well with the NLO charm quark distri- 
bution, as shown in Calderon's talk On the right-hand 
side of Fig. ^ we show the pr distributions at ^/S — 200 
GeV for the two sets of parameters in the total cross sec- 
tion curves on the left-hand side. The differences in the 
slopes are due to the different scales while the normaliza- 
tion difference is due to the choice of charm mass and the 
parton densities — the MRST densities generally give a 
larger cross section due to their larger as- However, the 
curves need to be scaled up by a factor of four to agree 
with the STAR normalization j4j , as may be expected from 
the total cross section results. 

To illustrate the effect of Peterson fragmentation, the 
Pt distribution with Ec = 0.06 in Eq. Q, is shown by the 
dot-dashed histogram in Fig. ^ The MRST parton den- 
sities are employed, along with m = 1.2 GeV, fi = 2m, 
as in the solid curve. At pt ~ 10 GeV, there is about a 
factor of 5 between the solid curve and the dot-dashed his- 
togram while there is a slight increase for pt < 2 GeV. If 
Cc is decreased by a factor of 10, making it more consistent 
with the Mellin space result the resulting dotted his- 
togram is rather similar to the charm quark distribution 



in the solid curve. A similarly hard fragmentation func- 
tion is used in the FONLL formalism, discussed in Frix- 
ione's talk [22], which matches fixed-order NLO terms to 
next-to-leading log, large pt resummation to produce an 
improved result for pT 3> m |25j . 

The shape of the charm quark pT distribution at = 
1.96 TeV also agrees quite well with the CDF data from 
the Tevatron |^ . Given the large discrepancy between the 
pQCD result and the STAR cross section, it might be sur- 
prising that the normalization is also in good agreement 
with the sum of the charged and neutral D data scaled 
to include Dg and Ac production. No total cross section is 
available because only charm hadrons with pT > 5 GeV 
have been measured so far. 

Factorization breaking for charm fragmentation has 
been suggested from studies of the xp distributions of 
e.g. and D~ production, particularly in Tr~ A inter- 
actions where the is leading relative to the _D+ since 
the shares a valence quark with the tt^ while the 
does not. Several mechanisms such as intrinsic charm jl8| 
and string drag have been proposed, both of which involve 
charm quark coalescence with spectators. Such comoving 
partons are naturally produced in a hard scattering. Al- 
though it is not intuitive to expect coalescence to work at 
high Pt, it seems to do so for charm. See the talks by Hwa 
|2n| and Rapp pi for discussions of coalescence models for 
light hadrons and charm. 

3 Hidden Charm 

We now turn to J/ip production in dA interactions at 
RHIC and the LHC. It is essential that the A depen- 
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dence be understood in cold nuclear matter to set a proper 
baseline for quarkonium suppression in A A collisions. The 
NA50 collaboration 1 , 27 has studied the J/^p A depen- 
dence and attributed its behavior to J/ break up by nu- 
cleons in the final state, referred to as nuclear absorption. 
However, it is also known that the parton distributions are 
modified in the nucleus relative to free protons. This mod- 
ification, referred to here as shadowing, is increasingly im- 
portant at higher energies, as we have recently emphasized 
|28| . In this section, we discuss the interplay of shadowing 
and absorption in d-f Au collisions at RHIC and in d+Pb 
collisions at the LHC. Previously, we calculated the effect 
of shadowing alone on the J/tp dA/pp ratio as a function 
of rapidity and impact parameter |2H|- The large cc to- 
tal cross section also has implications for the J/ip yield 
if J/ip's arise from cc recombination in a QGP |7|. Such 
a total cross section would suggest significant secondary 
J/ip production at RHIC, leading to enhancement rather 
than suppression in central collisions. There is no evidence 
for a strong regeneration effect in the PHENIX Au+Au 
data so far [SH]- 

Shadowing, the modification of the parton densities in 
the nucleus with respect to the free nucleon, may be taken 
into account by replacing /J in Eq. ^ by F^{x, n^,b, z) — 
PA{b, z)S^{A, X, fi^,b, z)fj{x, /i^). The density distribution 
of the deuteron is also included in these calculations but 
the small effects of shadowing in deuterium is ignored. 
We did not discuss the effect of shadowing on the charm 
Pt distributions because the effect at midrapidity is small 
and, on the logarithmic scale of the pt distributions, negli- 
gible. The J/tp is another story due to the PHENIX muon 
capability at forward and backward rapidity. As shown in 
Leitch's talk [30], although the PHENIX J/ip data are 
consistent with shadowing alone, the data are also consis- 
tent with nuclear shadowing plus a small absorption cross 
section of 1-3 mb, smaller than that currently obtained in 
SPS measurements |2I|. We have calculated J/ip produc- 
tion in the color evaporation model (CEM) using the same 
mass and scale as in cc production but cutting off the in- 
variant mass of the pair at 4m|). The calculations of the 
dA/pp ratios are done at LO to simplify the calculations. 
As shown in Fig. 13 the LO and NLO ratios are equivalent. 
We have now also implemented nucleon absorption in the 
calculation, showing the effect of several absorption and 
production mechanisms. 

To implement nuclear absorption on J/tp production in 
dA collisions, the dN production cross section is weighted 
by the survival probability, S^^^ 

S^'^^b, z) ^ cxp 1^- dz'pA{b,z')a,y,,{z' -z)^ (4) 

where z is the longitudinal production point and z' is the 
point at which the state is absorbed. The nucleon absorp- 
tion cross section, (Tabs, typically depends on where the 
state is produced in the medium and how far it travels 
through nuclear matter. If absorption alone is active, i.e. 

= 1, then an effective minimum bias A dependence 
is obtained after integrating S^^^ over the spatial coordi- 
nates. If S'^^'' = 1 also, adA = 2AapN- When S""^^ ^ 1, 



1 ' 




J 



Fig. 2. The J/tp pAu/pp ratio at 200 GeV. We compare 
the NLO (sohd histogram, MRST HO) and LO (soUd curve, 
MRST LO) results using m = /i/2 = 1.2 GeV with the EKS98 
parameterization. 

o-fjA = 2A"tTpAr where, if (Tabs is a constant, indepen- 
dent of the production mechanism for a nucleus of pA — 
PoO{Ra - b), a = 1- 9(Tabs/(167rr^), where rg = 1.2 fm. 
The contribution to the full A dependence in alxp) from 
absorption alone is only constant if (Tabs is constant and 
independent of the production mechanism The ob- 
served J/ip yield includes feed down from XcJ and ip' de- 
cays, giving 

S}Y^ib,z) = 0.58S}Y^^^,^{b,z) 

+ 0.35^;;'; (6, z) + O.US^'^^b, z) . (5) 

The J/ip may be produced as a color singlet, a color 
octet or in a combination of the two. In color singlet pro- 
duction, the final state absorption cross section depends 
on the size of the cc pair as it traverses the nucleus, allow- 
ing absorption to be effective only while the cross section 
is growing toward its asymptotic size inside the target. 
On the other hand, if the cc is only produced as a color 
octet, hadronization will occur only after the pair has tra- 
versed the target except at very backward rapidity. We 
have considered a constant octet cross section, as well as 
one that reverts to a color singlet at backward rapidities. 
For singlets, S'}^^.dir ^ ^t.j ^ -^V^'' "^^^^ octets, one 
assumes that Sf^^ ^j^, = S^f^ = 5'^V^ As can be seen 
in Fig. 13 the difference between the constant and grow- 
ing octet assumptions is quite small at large \/Snn with 
only a small singlet effect at y < — 2 and —5 at RHIC 
and the LHC respectively. Singlet absorption is also im- 
portant only at similar rapidities and is otherwise not dif- 
ferent from shadowing alone. Finally, we have also con- 
sidered a combination of octet and singlet absorption in 
the context of the NRQCD model, see Ref. |33 for more 
details. The combination of nonperturbative singlet and 
octet parameters changes the shape of the shadowing ra- 
tio slightly. Including the singlet contribution weakens the 
effective absorption. The results are shown integrated over 
impact parameter. The calculations use the EKS98 shad- 
owing parameterization since it gives good agreement with 
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Fig. 3. The J/t/; dA/pp ratios with EKS98 for VSj^ = 200 GeV (left-hand side) and for ^/SJ^ = 6.2 TeV (right-hand side) 
as a function of rapidity for (a) constant octet, (b) growing octet, (c) singlet, all calculated in the CEM and (d) NRQCD. 
For (a)-(c), the curves are no absorption (solid), CTabs = 1 (dashed), 3 (dot-dashed) and 5 mb (dotted). For (d), we show no 
absorption (solid), 1 mb octet/ 1 mb singlet (dashed), 3 mb octet/3 mb singlet (dot-dashed), and 5 mb octet/3 mb singlet 
(dotted). 



the trend of the PHENIX data. For resuhs with other 
shadowing parameterizations, see Ref. 32 . 

Several values of the asymptotic absorption cross sec- 
tion, (Tabs = I7 3 and 5 mb, corresponding to a = 0.98, 
0.95 and 0.92 respectively using Eqs. Q and (O are shown 
in Fig. O These values of (Tabs sue somewhat smaller than 
those obtained for the sharp sphere approximation. The 
diffuse surface of a real nucleus and the longer range of the 
density distribution result in a smaller value of (Tabs than 
a sharp sphere nucleus. There is good agreement with the 
trend of the preliminary PHENIX data ,33, for (Tabs = 0—3 
mb. We use a value of 3 mb in our further calculations to 
illustrate the relative importance of absorption and shad- 
owing. 

We now turn to the centrality dependence of J/tp pro- 
duction in dA collisions. In central collisions, inhomoge- 
neous shadowing is stronger than the homogeneous result. 
The stronger the homogeneous shadowing, the larger the 
inhomogeneity. In peripheral collisions, inhomogeneous ef- 
fects are weaker than the homogeneous results but some 
shadowing is still present. Shadowing persists in part be- 
cause the density in a heavy nucleus is large and approx- 
imately constant except close to the surface and because 
the deuteron wave function has a long tail. We also expect 
absorption to be stronger in central collisions. 

To study the centrality dependence of shadowing and 
absorption, we present the dA/pp ratios as a function of 
the number of binary NN collisions, A^coii, 

A^coii(6) = a'^^J d^sTAis)TB{\b~s\) 

where Ta and Tb are the nuclear thickness functions and 
crj^^ is the inelastic nucleon-nucleon cross section, 42 mb 
at 200 GeV and 76 mb at 6.2 TeV. In Fig.H we show the 
-^coii dependence for several representative rapidities, y = 



-2, and 2 for RHIC (left-hand side) and y -4, -2, 0, 
2 and 4 for the LHC (right-hand side). We have chosen an 
inhomogeneous shadowing parameterization proportional 
to the path length of the parton through the nucleus j2Hl • 
For results with other parameterizations, see Ref. |32| . 

The dependence of the RHIC ratios on -/Vcou is al- 
most linear, as seen on the left-hand side of Fig. 01 We 
do not show results for A^'coii < 1- The weakest A'cou de- 
pendence occurs in the antishadowing region, illustrated 
by the y = —2 result (dot-dashed curve). The trends of 
the ratios as a function of iVcou are consistent with the 
PHENIX data from the north muon arm {y = 2) and the 
electron arms {y = 0) but the preliminary PHENIX results 
from the south arm (y = —2) are much stronger than our 
predictions and, in fact, go the opposite way. The overall 
dependence on A'coU is stronger than that obtained from 
shadowing alone, described in Ref. 28 where inhomoge- 
neous shadowing effects depend strongly on the amount 
of homogeneous shadowing. Relatively large effects at low 
X are accompanied by the strongest b dependence. In the 
transition region around midrapidity at RHIC, the b de- 
pendence of the ratio dAu/pp due to shadowing is nearly 
negligible and almost all the iVcoii dependence at ?/ ~ can 
be attributed to absorption. The y = —2 results for color 
singlet production and absorption, in the antishadowing 
region, are fairly independent of A'^coU- 

On the right-hand side of Fig. 0] we present our in- 
homogeneous shadowing and absorption calculations for 
d+Pb collisions at ^/SNN = 6.2 TeV at the LHC. Results 
for y = ±4, in the range of the ALICE muon arm, are 
also included. Given that the rapidity range of the muon 
arm encompasses the crossover point where dPh/pp '-^ 1 at 
y —3.9, the centrality dependence of absorption alone 
could be determined and used to calibrate the inhomo- 
geneous shadowing effects. Note that is is only possible 
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coll coll 

Fig. 4. The dA/pp ratio as a function of A'coii for a growing octet with CTabs ~ 3 mb and the EKS98 parameterization. Left- 
hand side: results for y = —2 (dot-dashed), y — (dashed) and y = 2 (solid) at 200 GeV. Right-hand side: results for y = —4 
(dot-dot-dot-dashed), y = —2 (dotted), y — (dot-dashed), y = 2 (dashed) and y — 4 (solid) at 6.2 TeV. 



to reach y ~ —3.9 in ALICE by running Pb-|-d since the 
muon arm is only on one side of midrapidity. Both AL- 
ICE and CMS should be able to measure J/ip production 
at y = ±2 and 0. The results for y = A, ±2 and 0, all in 
the low X shadowing region, are rather closely grouped to- 
gether. This should not be surprising because the EKS98 
shadowing ratios shown on the right-hand side of Fig. |3 
are not very strong functions of rapidity. Note that the 
a;-axis scale is expanded relative to that of RHIC due to 
the larger crj^^ at 6.2 TeV. 



4 Summary 

In conclusion, the RHIC d+Au data on open charm and 
J/ip are beginning to come into their own. While the 
QCD calculations agree well with the shape of the STAR 
Pt distributions, they underestimate the reported total 
cross section. In contrast, the J/tp cross section is in rel- 
atively good agreement with QCD predictions and the 
agreement of the minimum bias data with calculations in- 
cluding shadowing and nucleon absorption is quite good. 
The agreement of the J/ip calculations with the prelimi- 
nary PHENIX data is generally quite good, except for the 
dependence of the y ~ —2 results on iVcoii. 
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